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ABSTRACT 

We present the results of our survey toward the protoplanetary nebula CRL 
618 for several large, highly saturated, oxygen bearing organic molecules of bio- 
logical importance including acetaldehyde (CH 3 CHO), acetic acid (CH3COOH), 
dimethyl ether (CH3OCH3), ethanol (CH 3 CH 2 OH), formic acid (HCOOH) and 
methyl formate (HCOOCH3); large carbon chain molecules including methyl 
cyanide (CH 3 CN), methylcyanoacetylene (CH3C3N), cyanoacetylene (HC 3 N), 
cyanodiacetylene (HC 5 N), and C 6 H; and finally smaller molecules including 34 SO, 

S 0 2 , OC^S and MgNC. No biologically important organic molecules were de- 
tected. However, we report the first interferometric detections of CH 3 CN and 
vibrationaUy excited HC 3 N and HC 5 N toward this source. The temperature and 
distribution of CH 3 CN toward CRL 618 indicates it is formed in the outer enve- 
lope surrounding the UC HII region. Furthermore, the P-Cygni line profile and 
corresponding channel maps of vibrationaUy excited HC5N supports its distri- 
bution in the extended envelope expanding radially from the central star. The 
detection of vibrationaUy excited HC3N confirmed the temperature structure and 
column density of HC3N in the inner envelope found by Wyrowski and coUeagues 
( 2003 ). Finally, our observations clearly indicate that CRL 618 is a good source 
of large carbon chain species but is a very poor source to detect or produce 
organic species of biological importance. 
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1. Introduction 

Recently, there have been several surveys for complex molecular species toward mas- 
sive star-forming regions including W75N, G45.47+0.05, G19.61-0.23 (Remijan et al. 2004a) 
and Sgr B2(N-LMH) (Hollis et al. 2004; Friedel et al. 2004, Snyder et al. 2002). These 
massive star-forming regions all have very similar physical conditions. First, they are rich 
in interstellar dust. Many molecular species, even H 2 , appear to have their origins on the 
icy surfaces of dust grains (Hollenbach & Salpeter 1971). Dust becomes more important 
as more complex molecules are found in the interstellar medium (ISM). The abundances 
of large molecules in several regions are much larger than the predicted abundances from 
ion-molecule gas phase chemistry' models (Remijan et al. 2002, Mehringer & Snyder 1996). 
Thus, grain surface chemistry has been invoked to account for the large abundance of com- 
plex interstellar molecules. Second, the high density of gas and dust limit the amount of UV 
flux. This is important because UV photons destroy complex molecules (Ehrenfreund et al. 
2001). Finally, the hot core regions have high temperatures. It is generally believed that 
molecules formed via grain surface reactions in a cold collapse phase are then evaporated by 
thermal heating and thereafter deposited into the gas phase. Typically, if the temperature 
of dust grains exceed 100 K, the mantles evaporate and species that were frozen on the grain 
can now escape into the gas phase (Tielens & Hagen 1982; Hasegawa, Herbst & Leung 1992). 
Furthermore, an interesting warm gas phase chemistry may take place after the molecules 
are evaporated leading to the production of even more complex molecules. Thus, these re- 
gions appear to be ideal environments for complex interstellar chemisty to take place. It is 
logical to assume that similar regions rich in dust, high density, and near a source of high 
temperature, i.e. O and B stars or ultracompact (UC) HII regions, should also contain the 
necessary chemistry to synthesize large molecular species. 

Another astronomical environment that appears to fit the above selection criteria are the 
extended envelopes of protoplanetary nebulae (PPNe). The PPN phase of stellar evolution 
occurs at the end of the asymptotic giant branch (AGB) stage, where substantial mass loss 
is taking place, and before the formation of the white dwarf, after which the temperature 
becomes high enough to ionize the extended envelope and start the formation of a planetary 
nebula (PN). This phase of stellar evolution is typically very short, lasting only about ten 
thousand years. However, during that time a very interesting chemistry may take place in 
the extended envelope of the post- AGB star. Thus, the extended envelopes of PPNe may be 
ideal chemical factories to produce large molecules and molecules of biological importance. 

CRL 618 (AFGL 618, IRAS 04395+3601) is a carbon-rich, evolved PPN with a thick 
molecular envelope surrounding a B0 star and an UC HII region (Cernicharo et al. 2001a). 
However, the complete physical structure of CRL 618 is very complex. CRL 618 has been 
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observed in several molecular species at different resolutions and wavelengths. Some of its 
distinct structures include optical high velocity bipolar outflows (Trammell 2000), a cen- 
tral UC HII region, a low velocity expanding torus of molecular emission and an extended 
expanding AGB envelope (Sanchez Contreras & Sahai 2004). In addition to several indi- 
vidual structures, the CSE itself also shows a complex structure and is predicted to exhibit 
a complex distribution of several different molecular species (Woods et al. 2003). Thus, by 
observing the molecular envelope in several different molecular species at different excita- 
tion temperatures, we can obtain a better understanding of the temperature and density 
stucture of the CSE. Furthermore, the interactions of the high and low velocity outflows 
from the embedded evolved star may also contribute significantly to the chemistry in the 
CSE. For example, within the high velocity bipolar outflow, Ueta, Fong & Meixner (2001) 
discovered distinct near-IR “microstructures”. They note that the positions of these “mi- 
crostructures” occur in shocked regions where the molecular clumps are interfacing with the 
ambient circumstellar shell. Shocked regions also often trace areas of interesting interstellar 
chemistry that lead to the formation of complex molecular species. For example, Liu et al. 
(2002) detected formic acid (HCOOH) emission from shocked gas toward the Orion Com- 
pact Ridge. Also, Hollis et al. (2004), detected emission from propenal (CH 2 CHCHO) and 
propanal (CH 3 CH 2 CHO) toward Sgr B2(N). The formation and widespread distribution of 
these species are most likely due to shocks in the Sgr B2(N) region that provided the en- 
ergy necessary for the endothermic reactions to form these species and eventually distribute 
them into the gas phase. Thus, it appeared that CRL 618 contains the physical conditions 
necessary to produce complex molecules. 

Two recent surveys were performed toward CRL 618 utilizing the Infrared Space Ob- 
servatory (ISO) and the IRAM 30 m radio telescope (Cernicharo et al. 2001a, b). These 
observations made the first detections in the circumstellar medium of the polyacetylenic 
chains C 4 H 2 and C 6 H 2 and made the first detection of circumstellar benzene (C 6 H 6 ). They 
also detected the infrared bands of NH 3 and C 2 H 4 and made the first detection in the cir- 
cumstellar medium of the methylpolyyenes CH 3 C 2 H and CH 3 C 4 H. Obviously, there is a rich 
chemistry available to form large carbon-chain molecular species, including HC 3 N and HC5N 
(Fukasaku et al. 1994). Bujarrabal et al. (1988) also compared the column densities of sev- 
eral molecular species measured toward CRL 618 to what was measured toward other PPN, 
specifically IRC+ 10216 and CRL 2688, as well as interstellar molecular clouds, specifically 
Orion A and TMC-1. What they found was that the abundances of many molecular species, 
including HC 3 N and HC5N, were more similar to the abundances found in the ISM rather 
than the other PPNe. Furthermore, it was the first time HCO + was detected in a carbon- 
rich object. Thus, they suggested that CRL 618 was at a particular point in its evolution 
where the chemistry occuring in the CSE was more likely similar to the chemistry present 
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in the ISM than other, more evolved, PPNe. There have been detections of small oxygen 
bearing molecules toward CRL 618 including CO, H 2 0 and H 2 CO (Cernicharo et al. 1989). 
However, there has not been a detection of any large oxygen bearing organic molecule toward 
this source even though CRL 618 appears to be the ideal CSE environment to search for 
oxygen bearing organic molecules as well as other large molecular species. 

Thus, to investigate whether the chemistry present toward CRL 618 can produce large 
oxygen bearing organic molecules, we used the Berkeley-DIinois-Maryland Association (BIMA) 
Array 5 to carry out a large molecule survey toward this region. Included in the survey were 
transitions of formic acid (HCOOH), acetic acid (CH 3 COOH), methyl formate (HCOOCH3), 
acetaldehyde (CH3CHO) as well as methyl cyanide (CH 3 CN), cyanoacetylene (HC3N) and 
cyanodiacetylene (HC5N). The observations of these large molecules not only helped con- 
strain the chemistry present toward CRL 618 but also the temperature and column density 
of large molecular species in the surrounding circumstellar environment. 


2. Observations 

The large molecule survey took place from 2002 October to 2003 February utilizing the 
B and C configurations of the BIMA Array. The minimum baseline for these observations 
was ~8 m, and the maximum baseline was ~150 m, providing angular resolutions of ~2"- 
13". Table 1 lists the observational parameters of our survey toward CRL 618. Column (1) 
lists the observation date. Column (2) lists the source used to calibrate the antenna based 
gains and column (3) lists the flux density calibrator. The absolute amplitude calibration 
of these sources is accurate to within ~20%. We find a flux density of the source of 1.86 
Jy at the continuum emission peak of our observations at 108 GHz which is very similar 
to the flux density measured by Wyrowski et al. (2003) of 1.8 Jy measured with the IRAM 
30 m radio telescope at 109.2 GHz. Column (4) lists the observing frequency of the upper 
sideband and column (5) lists the array configuration of the observation. Finally, column 
(6) lists the corresponding beamsize. In order to obtain the best signal-to-noise at each 
observing frequency, all the data from the source with multiple tracks and in multiple arrays 
were inverted in u-v space. The resulting synthesized beamsize (5".8x5".l near 108 GHz) 
was determined by the region of densest u-v coverage and the 1 a rms noise level of the 
final maps was determined by selecting a region free from line and continuum emission. The 


5 Operated by the University of California, Berkeley, the University of Illinois, and the University of 
Maryland with support from the National Science Foundation. 
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passbands were calibrated automatically during data acquisition 6 and this method has been 
quite satisfactory and has not generated spurious features. 

Table 2 lists the molecular parameters of the 44 transitions of the 15 molecular species 
searched for in our survey. Column (1) lists the species; Column (2) lists the transition. 
Column (3) lists the calculated rest frequency and its 2 o standard deviation. Column 
(4) lists the product of the line strength and the square of the relevant dipole moment; 
and column (5) lists the upper state energy level of each transition. Finally, the estimated 
rotational partition functions for each molecular species are given in column (6). The spectral 
windows containing these transitions had a bandwidth of 50 MHz and were divided into 128 
channels for a spectral resolution of 0.39 MHz per channel. All data were combined, imaged, 
and self-calibrated using the MIRIAD software package (Sault, Teuben & Wright 1995). The 
1 a rms noise level of the final maps was used to set upper limits to the column densities of 
the undetected molecular species searched for in our survey (§4.4). 


3. Results 

Figure 1 shows the four transitions of vibrationally excited HC 5 N located in the pass- 
bands of our survey. The HC5N transitions and corresponding rest frequencies were taken 
from a new survey by Yamada et al. (2004) and from the Cologne Database for Molecular 
Spectroscopy (Muller et al. 2001). In Figure la) we show the first interferometric detection 
of vibrationally excited (u u =l) HC 5 N toward CRL 618. The spectral line labels located in 
the upper left of each panel correspond to the rest frequency of each HC 5 N transition for 
a v L sr of -27.5 km s -1 (dashed line). The 1 a noise level is seen at the left of each panel. 
The classic P-Cygni line profile indicates that HC5N is present in an expanding envelope of 
material. Also included in Figure la) are the spectral line labels of C 6 H and ^SO. These 
features will be discussed in more detail (§§4.4 & 4.5). Figures lb)-ld) show the spectra of 
higher vibrational transitions («u = 3 and 4) that were not detected beyond our 1 a detec- 
tion limit. From these non-detections, we were able to set a upper limit on the temperature 
of this species in this region. Finally, in Figure Id), we label the location of MgNC. The 
non-detection of this feature will also be discussed in more detail (§4.4). 

Figure 2 shows the spectra of the CH3CN and vibrationally excited HC3N detected in 
our survey. Figure 2a) shows the first detection of the J = 6 — 5 rotational transitions of 
CH 3 CN. All the CH 3 CN transitions are seen in emission. Between the K = 2 and K = 3 
lines of CH 3 CN lies the J = 12 — 11, (0,1,0,3) 1+ vibrational transition of HC 3 N seen in 


6 A technical description of this can be found at http://astron.berkeley.edu/~plambeck/technical.html 
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absorption. The measured v LS r of the absorption trough is -27.5 km s -1 which is consistent 
with previous detections of vibrationally excited HC3N (Wyrowski et al. 2003). Finally, 
Figure 2b) shows the spectrum of 3 vibrationally excited transitions of HC 3 N. Once again, 
each of these transitions lie at a vlsr of -27.5 km s -1 . The spectral line labels located in 
the upper left of each panel correspond to the rest frequency of the HC 3 N J = 12 — 11, 
(0,1,0,3) 1+ , and the HC 3 N J = 12 — 11, (0,0, 0,4)° transitions 7 , respectively. Table 3 lists the 
molecular species that were detected toward CRL 618 at or near the systemic velocity (-27.5 
km s -1 ). Least square Gaussian fits were made for each spectral line in order to obtain 
the radial velocities and integrated fine intensities for the detected transitions. In Table 
3, column (1) lists the calculated rest frequencies of the identified species. Columns (2) 
and (3) list the molecular species and associated transitions, respectively. Column (4) lists 
the v L sr corresponding to the emission or absorption components of each spectral feature. 
Finally, column (5) lists the fitted integrated line intensity of each detected transition and 
the associated 1 a uncertainty. 

Figure 3 shows the HC5N velocity-channel maps toward CRL 618. The spectral line 
channel maps displayed are a combination of the B and C array data. The v LS r intervals are 
indicated in km s -1 units at the top right of each panel. The synthesized beam of 5."8x5."l 
is shown at the bottom left of each velocity-channel map. The lowest contour plotted is 
0.08 Jy beam -1 which corresponds to 2 a. The absorption peaks near a vlsr of -34 km 
s -1 and the emission peaks near -24 km s -1 . Both the absorption and emission peaks are 
coincident with the high resolution B-array 3 mm continuum shown in the last panel of the 
figure (synthesized beam of 3."0xl."9). Thus, it is apparent from both the P-Cygni profile 
and velocity-channel maps that HC5N is contained in an expanding envelope and not in the 
bipolar outflows. 

Figure 4 shows the HC 3 N velocity-channel maps toward CRL 618. As in Figure 3, the 
data displayed are a combination of the B and C array data and the vlsr intervals are 
indicated at the top right of each panel. In this case, we see no emission features of HC 3 N. 
This is not unexpected because all of the high energy, high vibrational state transitions 
(i.e. u 5 = 1; v 7 = 3) detected by Wyrowski et al. (2003) were seen primarily in absorption. 
The transitions detected in our survey were a higher vibrational state (v 7 =3 & 4) than 
the transitions of the Wyrowski et al. (2003) survey but the upper state energy levels were 
comparable to some of the higher energy transitions detected in their survey. In contrast, 
the ground state and first vibrational states in the Wyrowski et al. (2003) survey were 
seen primarily in emission. Once again, we see the absorption is coincident with the 3 mm 


7 The notation given for the vibrational excited transitions of HC 3 N is defined as (v4,v 5 , i/ 6 , V7) 1 , where v 
is the vibrational quantum number and l is the l - type doubling quantum number (Thorwirth et al. 2003). 
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continuum and presumably all the vibrationally excited HC 3 N is contained in the CSE. 

Finally, we note that we did not detect any emission or absorption features from any of 
the other molecular species listed in Table 2 beyond our 1 a detection limit. These include 
all the oxygen bearing organic molecules as well as CH 3 C 3 N, C 6 H, and all of the smaller 
molecular species. However, it does appear that we are starting to detect the N = 9-8 
transition of MgNC in Figure Id). Higher signal-to-noise interferometric observations are 
necessary to confirm this detection. Based on the 1 a detection limit, we were able to place 
upper limits on the species that were not detected. All column densities and upper limits 
are listed in Table 4 (§4.4). In §4.5 we discuss these upper limits along with the predicted 
fractional abundances based on the gas phase formation models of Woods et al. (2003) and 
of the UMIST database (Le Teuff et al. 1999). 


4. Analysis and Discussion 

One of the most important physical parameters to determine in either an interstellar 
or circumstellar environment is the total beam averaged column density of a molecular 
species. There are three primary procedures used in calculating the total beam averaged 
column density of a particular molecular species, each of which are described in detail in 
Remijan et al. (2004b). For our purposes, we will utilize the rotational temperature diagram 
method (RTDM) for the CH 3 CN and HC5N emission, and to assume local thermodynamic 
equilibrium (LTE) for all other molecular species. From Remijan et al. (2003), the total 
beam averaged column density is given as: 


N t = 2.04 x C T 


/A/du f /(r e ) 

e a e b i(T e ) - i b \ v 3 < StjiJt 2 > 


Q e (B u /T e ) 

Wrv x 10 20 cm" 2 = C T Np in cm -2 . (1) 


In equation (1), N T is the total column density for any optical depth and Nft tn is the total 
column density for the optically thin case. C T is the optical depth correction factor from 
Goldsmith & Langer (1999). / A Idv (Jy beam -1 km s _1 ) is the integrated line intensity 
given in Table 3; 9 a and 0 b are the FWHM beam sizes (arcsec), I(T e ) is the specific intensity 
of the source at T e (Jy beam -1 ), and I b is the specific intensity of the total background 
continuum contribution (Jy beam -1 ). Q rv is a combination of the rotational and vibrational 
partition functions. The vibrational part of the partition function needed to be applied 
to both the HC5N (Muller et al. 2001) the HC 3 N (Wyrowski et al. 2003) column density 
calculations. E u is the upper state energy level (K), T e is the rotational temperature (K), 
v is the rest frequency (GHz), and S^^i 2 is the product of the total torsional-rotational line 
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strength and the square of the electric dipole moment (Debye 2 ). 


4.1. Methyl Cyanide (CH 3 CN) 

Methyl cyanide (CH 3 CN) has been used as an excellent tracer of the physical and 
chemical conditions toward high mass star forming regions of the interstellar medium (ISM) 
(Remijan et al. 2004b; Sutton et al. 1986). CH 3 CN is a symmetric top molecule which has 
properties that makes it an ideal probe of kinetic temperature and column density. In a 
symmetric top, the rotational levels are labeled by J, the total angular momentum quantum 
number, and K , the quantum number for the projection of angular momentum along the 
symmetry axis (Cummins et al. 1983). By selection rules, radiative transitions do not change 
the value of K. However, collisional transitions are permitted across K ladders. Therefore, 
the population of one K ladder relative to another is largely determined by collisions and 
hence, is closely related to the kinetic temperature and density (Solomon et al. 1971). The 
first reported detection of CH 3 CN toward CRL 618 was by Bujarrabal et al. (1988) using 
the IRAM 30m radio telescope. Also using the IRAM 30m, Guelin & Cemicharo (1991) 
used several millimeter transitions of CH 3 CN to examine the physical conditions toward 
IRC+10216. From their observations, they found a rotation temperature of T ro t = 16 K and 
a column density of 6 xl 0 12 cm -2 toward IRC+10216 (Guelin & Cemicharo 1991). Figure 5 
shows the results of the RTDM toward CRL 618 using the J = 6 — 5 transitions. Because the 
K = 0 & 1 transitions are severely blended, we could not resolve the lines into two separate 
Gaussians, thus we included the pair as one point on the rotational temperature diagram by 
separating the intensity by their appropriate statistical weights and then averaging. Using 
the procedures outlined in Remijan et al. (2004b) we found a best fit rotational temperature 
of 88(44) K and a total beam averaged column density of 3.2(3) xlO 14 cm -2 . Assuming an 
H 2 column density of ~1.3xl0 23 cm -2 (Woods et al. 2003), we find a fractional abundance 
of Ach 3 cn=2.4(3)x10 -9 . This value is very close to the predicted peak value of ~ 10 -8 from 
the Woods et al. (2003) model at a radial distance of ~9.1xl0 15 cm from the exciting star. 
Thus, it is clear from our observations that CH 3 CN appears to reside in the outer regions 
of the expanding envelope. Furthermore, transitions of vibrationally excited CH 3 CN v% = 1 
were below our 1 a detection limit, which again may imply that CH 3 CN exists in the outer 
envelope. 
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4.2. Cyanoacetylene (HC 3 N) 

Cyanoacetylene (HC3N) is a carbon chain molecule that has been extensively studied 
toward CRL 618 in both the ground and excited vibrational states (see e.g. Pardo et al. 
2004; Wyrowski et al. 2003 and references therein). So, from our data, we were able to check 
the previously measured physical conditions and help further constrain the physical and 
chemical models of CRL 618. In Figure 2), we show the detections of 4 vibrationally excited 
transitions of HC 3 N. Also shown in these spectra are the best fit modelled spectra to the 
data given as a solid line. In order to calculate the best fit spectra, we used the spectral line 
parameters given in Table 2 and used a FWHM line width of Av = 5 km s -1 based on the fits 
in Table 3. Models were constructed using the temperature and upper state column density 
as input parameters. We varied the temperature between 400 and 600 K and the upper state 
column density between 10 16 and 10 18 cm -2 . We then calculated the x 2 between the model 
and the data and found the minimum x 2 - The best fit models (solid lines) are shown in 
Figure 2) for a temperature of 580(150) K and a total HC3N beam averaged column density 
of 4.7(6) xlO 17 cm -2 . The errors given for the temperature and total column density are 3 
a errors based on the fit, thus giving a 99% confidence interval. The only line that does not 
seem to fit the model very well is the J=12-ll (0,0, 0,4) 4 transition at 110.555 GHz (Table 2). 
This line may be blended with the 15 4 ^-14 5i2 transition of C-C3H2 at 110.558(9) GHz (Pickett 
et al. 1998) or possibly from another unknown circumstellar species. The values of T rot and 
Nr are very similar to the values determined by Wyrowski et al. (2003) of T=520(80) K and 
a total HC 3 N column density between (3-6) xlO 17 cm -2 . In order to find the total HC 3 N 
column density, we needed to include the vibrational part of the partition function. Using 
Figure 9) in Wyrowski et al. (2003), for a temperature of 580 K, we find Q rt ,=1.8xl0 5 . 
Also, in equation (1) we needed to account for the optical depth and the intensity of the 
background continuum. For the HC 3 N lines, we find r=ln(T c /T c -Tx,) = 0.25 which gives an 
optical depth correction factor of C T = 1.13 and using the continuum temperature and the 
excitation temperature of 580 K, the column density correction due to /(r ei )/(/(T aE ) — h) 
is ~1.06 which gives a total correction of ~1.20 for a temperature of 580 K. Assuming an 
H2 column density of ~1.3xl0 23 cm -2 (Woods et al. 2003), we find a fractional abundance 
of Ahc 3 h=3.6(5) x 10 -6 , which match the models of Woods et al. (2003). The average value 
of the HC 3 N fractional abundance in the envelope predicted by the model is ~10 -6 and the 
model predicts a peak value of 1.6xl0 -6 at a radial distance of ~10 16 cm from the central 
star. Thus, it appears that the models are predicting the relative abundances of this species 
very well toward CRL 618. 
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4.3. Cyanodiacetylene (HC 5 N) 

The first reported detection of HC 5 N in a PPN was by Winnewisser & Walmsley (1978) 
toward IRC+10216. They observed the J — 9 — 8 rotational transition near 24 GHz. From 
their observations, they found the HC 5 N emission to be quite extended over their 40" beam 
and measured a HC5N column density of 4xl0 14 cm -2 . However, they did not detect any 
emission from the un=l vibrational mode. These observations were later followed up by 
Jennings & Fox (1982) who determined an HC5N rotational temperature of 13^4 2 K, which 
was in agreement with the determined HC 3 N rotational temperature of 12 K measured at the 
time. The first detection of HC 5 N in CRL 618 was by Bujarrabal et al. (1988). In this case, 
they assumed the HC 5 N emission to be much smaller than their beam of ~15". From that, 
they found an HC 5 N rotational temperature of ~81 K. Since then, several other transitions 
of HC 5 N have been observed toward other PPNe (see e.g. Fukusaka et al. 1994). The first 
detection of vibrationally excited transitions of HC 5 N reported in CRL 618 was by Thorwirth 
(2001). Later observations followed by Cernicharo et al. (2001a) and Wyrowski et al. (2003). 
However, ours are the first data of vibrationally excited HC5N taken with an interferometer 
and the first maps of vibrationally excited HC 5 N toward CRL 618. In Figure 1, we showed 
the first interferometric detection of the J = 40 — 39, un=l transition of HC5N. No higher 
vibrational states (vn=3 or 4) were detected (Figures lb-d)). However, from the 1 a noise 
levels, we were able to plot the upper limits on a rotational temperature diagram of HC 5 N. 
Plotting these data on a rotational temperature diagram, we were able to determine an upper 
limit to the HC 5 N rotational temperature of T rot = 114 K. Using the measured vibrationally 
excited partition function from Muller et al. (2001) for HC5N, we were able to calculate the 
total rotational- vibrational HC 5 N partition function for 114 K. We find Q r „ = ~950. From 
the rotational temperature diagram, we find a total beam averaged column density of HC5N 
of ~10 15 cm -2 and a fractional abundance of Xhc 5 n=7.6x 10 -9 . This column density is about 
an order of magitude lower than the predicted value of 1.7xl0 16 cm -2 from the Woods et al. 
(2003) model. However, the column density is in the range of measured values toward t h i s 
source from previous observations (Fukasaku et al. 1994, Cernicharo et al. 2001a). 

Also, looking at the velocity-channel maps in Figure 2, we see the HC 5 N emission is 
coincident with the 3 mm continuum and that no HC 5 N emission is seen in the bipolar 
outflow above our 3 a detection limit cutoff. The classic P-Cygni line profile indicates 
a spherically symmetric expansion with an expansion velocity of ~5-7 km s -1 . The slow 
expansion velocity would occur at the outer edges of the envelope as the expanding material 
begins to interact with the remnant AGB wind. Also, the expansion velocity of HC 5 N is 
close to the measured expansion velocity of HC3N (~10 km s -1 ) found by Pardo et al. (2004) 
and it matches the expansion velocity used by Woods et al. (2003) in their formation models. 
Thus from the spectral line profile and mapping data, it appears that HC 5 N is found in the 
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extended circumstellar envelope. 

To summarize, all the detected species are found in the extended envelope. In the 
case of CH 3 CN, the observed transitions are low in upper state energy and optically thin. 
Thus, we are probing the gas in the very extended envelope and should only see emission 
features. On the other hand, the HC3N transitions have high upper state energies. The 
gas sampled in this case must lie deeper in the envelope where the column density is high 
which leads to a high opacity and the lines are seen primarily in absorption. Transitions 
of HC3N at a lower energy state with respect to the transitions we observed are seen by 
other authors in emission and as P-Cygni profiles. This implies these authors are sampling 
the outer envelope where the gas is at a lower temperature but the column density is high 
enough where there are competing effects of emission and absorption (see e.g. Pardo et al. 
2004). Finally, for HC5N, the temperature implies the gas resides in the outer envelope, 
intermediate between the vibrational HC 3 N and lower energy CH3CN transitions, and the 
column density is high enough where, once again, there are competing effects of emission 
and absorption. The hypothesis based on the formation models of each species (Woods et al. 
2003) is that CH 3 CN, HC3N and HC5N are probably contained throughout the CSE of CRL 
618. In order to fully investigate that assumption, observations of transitions that sample a 
wide range of upper state energy levels of each species are necessary. 


4.4. Column Density Upper Limits of the Remaining Molecular Species. 

Of the remaining species listed in Table 1), we only detected emission from ^SO (Figure 
la). Because this line was near the edge of the spectral window, we were not able to make a 
Gaussian fit or a calculation of the total ^SO beam averaged column density. However, this 
is the first detection of this isotopic species of SO in a PPN. Further observations of 34 SO 
and 32 SO will be needed in order to determine the isotopic ratio and compare the column 
densities with other detected O bearing species including CO and H2O toward CRL 618 
(Herpin & Cemicharo 2000). 

We also note that we did not reliably detect any emission from the TV = 9—8 transition of 
magnesium isocyanide (MgNC). This species was recently detected using the Arizona Radio 
Observatory (ASO) 12 m and the IRAM 30 m telescopes via five rotational transitions at 
95.4, 131.2, 143.2, 151.1, and 167.0 GHz. The column density of MgNC from that survey 
was found from a rotational temperature diagram to be 2.4xl0 12 cm -2 with a T rot = 21 K 
(Highberger & Zuirys 2003). They note the derived rotational temperature is consistent 
with the kinetic temperature of ~50 K (Herpin et al. 2002). Using this column density, we 
can predict the expected line intensity for our observations. Using a rotational temperature 
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of 21 K, we predict an intensity of ~0.035 Jy beam -1 for our BIMA observations. If the 
temperature is closer to the CH 3 CN temperature of 88 K, the predicted intensity drops 
to ~0.022 Jy beam -1 . The average rms noise of the combined B and C array dataset 
containing the MgNC transition is ~0.04 Jy beam -1 . Thus, even with the highest predicted 
line intensity, we would have only detected this transition at the ~1 a level. In Figure Id), 
we labeled the position of the MgNC transition and it appears the line is beginning to be 
resolved out of the noise. However, higher signal-to-noise interferometric observations are 
necessary to claim this as a detection. 

No other molecular species was detected toward CRL 618 beyond the 1 a detection 
limit. However, we were able to set upper limits to the column densities of those species 
based on the 1 a rms noise level (~0.04 Jy beam -1 ) and an average line width of ~5 km 
s -1 and assuming a temperature of 88 K found from the CH 3 CN data. Table 4 summarizes 
the column densities of all molecular species including upper limits. Column (1) lists the 
molecular species observed toward CRL 618. Column (2) lists the frequency of the transition 
used to calculate the upper limit to the total column density (MHz). Using equation (1), 
an average beam of 25 arcsec 2 , and the spectroscopic parameters and partition functions in 
Table 2, we were able to calculate the total beam averaged column densities of each molecular 
species. Column (3) lists the total beam averaged column density (cm -2 ). Finally, Column 
(4) lists the fractional abundance of each molecular species with respect to H2 assuming an 
H2 column density of 1.3X10 23 cm -2 . 


4.5. Formation and Detection of Oxygen Bearing Organic Biomolecules 

in Protoplanetary Nebulae 

Our detections of CH 3 CN, HC 3 N and HC5N strongly suggest that the chemistry present 
toward CRL 618 is conducive to producing large carbon-chain molecules and the formation 
models that predict the abundance of these species also match the measured abundances 
very well. The only large carbon-chain molecule that was predicted to be present in a large 
abundance toward CRL 618 but was not detected was CgH. Woods et al. (2003) predict 
a total column density of 2.5 xlO 15 cm -2 and our measured upper limit was two orders of 
magnitude lower at 2.5 x 10 13 cm -2 . Also, there is no formation model that takes into account 
methylcyanoacetylene (CH 3 C 3 N) and our measured upper limit is 3.3 xlO 13 cm -2 . However, 
the models of Woods et al. (2003) predict the column densities of CH 3 C,4H and CeH 5 CN 
to be 1.5 xlO 15 cm -2 , which are presumed to be more complex species than either C 6 H or 
CH 3 C 3 N. Thus, if the formation models are correct, we should have easily detected C 6 H. If 
the predicted column density of CH 3 C 3 N is similar to either CH 3 C4H or C 6 H 5 CN, we should 
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have also easily detected CH3C3N. The fact that neither species was detected beyond our 1 a 
detection limit implies the calculated upper limit of ~10 13 cm -2 for both C 6 H and CH3C3N 
must be taken into account in future formation models. 

The non-detections of very favorable transitions of oxygen bearing organic molecules 
such as HCOOH, HCOOCH3 and CH3COOH, however, imply that the chemistry of PPNe 
may not be conducive to form large organic molecules of biological interest (biomolecules). 
In fact, the only biomolecules that have been detected toward a PPN are formaldehyde 
(H 2 CO) (Cemicharo et al. 1989) toward CRL 618 and a tentative but unverified detection 
of methanol (CH3OH) toward IRC+10216 (Latter &; Chamley 1996). There are numerous 
chemical models that have been utilized to try and explain the formation and measured 
abundances of oxygen bearing organic molecules in the ISM and CSE. These models include 
many different types of energetics including shocks and thermal heating to drive the predicted 
chemistry. For our purposes, we will only include a brief overview of some models and how 
they may apply to CSEs. 

In the case of forming CH3OH and H 2 CO in a CSE, shock chemistry is invoked to 
explain the formation of H 2 CO and CH3OH. It is believed that both species would form 
in the post-shock gas after CO dissociates. For CH 3 OH, CH4 + OH or CH 3 4- H 2 0 could 
form CH3OH in the postshock gas (Latter & Chamley 1996; Hartquist et al. 1995). H 2 CO 
would form in a similar manner with the reaction of O 4- CH 3 and then the release of a 
H atom (Cemicharo 2004). Another possible formation mechanism is H addition after the 
dissociation of CO in the post-shock gas. 

Grain-surface reactions are also proposed to account for the abundance of oxygen bear- 
ing molecules. For example, Watanabe et al. (2004, and references therein) have shown that 
H 2 CO and CH3OH can form from both CO and C0-H 2 0 ice sprayed with cold H. From this 
experiment they concluded that for an ISM environment with 1 H atom cm -3 and a temper- 
ature of 10 K, a timescale of 10 6 yr would be necessary to match observed CO/CH3OH and 
H 2 CO/CH 3 OH ratios in molecular clouds. Subsequent H addition has also been theorized 
to produce complex oxygen bearing organic molecules including propenal (CH 2 CHCHO) 
and propanal (CH3CH 2 CHO) (Hollis et al. 2004). To date, however, CH 3 OH has not been 
detected toward CRL 618 and the timescale of the PPN phase is only about ten thousand 
years. It is apparent that this formation mechanism would be too slow to produce detectable 
amounts of CH 3 OH in CRL 618. 

Finally, there are also gas phase reaction schemes that attempt to predict the abundance 
of oxygen bearing organic molecules. In order to investigate whether gas phase reactions may 
be responsible for the formation of oxygen bearing organic molecules in the CSE of CRL 
618, we utilized the gas phase reactions of the UMIST database (Le Teuff et al. 1999) and 
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compared the predicted fractional abundance to the upper limits measured from our survey. 
The limitation that we found to this procedure is the limited number of reactions to form the 
larger oxygen bearing organic molecules. However, we can discuss two molecular formation 
models and how they compare to our measured upper limits. Furthermore, as a constraint to 
our analysis, we only included formation models that contained molecules already detected 
toward CRL 618 or molecules that were structurally similar to other species but have not 
yet been detected. 

We adopt a measured UV flux of ~3xl0 6 photons s _1 , a nn 2 abundance of 10 7 cm -3 
and a temperature of 100 K (Herpin & Cemicharo 2000). The proposed gas phase reaction 
formation of HCOOH is HCO + + H2O— ►HCOOHj and the destruction mechanism HCOOH 
4- hz/— >HCO + OH. HCO + and H2O have been detected in the CSE of CRL 618 and their 
fractional abundances have also been measured. The calculated fractional abundance of 
HCOOH based on the formation model is ~5xlQ -11 . The calculated upper limit based on 
our observations is Xhcooh~2x10 -9 . Thus, the combination of a lower abundance of H2O 
and HCO + and a higher UV ionization rate suggests that doubly-oxygenated molecules have 
abundances that are still too low to detect. No other oxygen bearing organic molecule we 
searched for had any formation mechanism based on species already found in the CSE of 
CRL 618 or formation models listed in the UMIST database. 

However the situation is less clear for the singly oxygenated species such as CH3CHO. 
The proposed gas phase reaction formation of CH 3 CHO is O + C2H5-4CH3CHO -I- H and 
the destruction mechanism CH3CHO + hv — ►CO + CH 4 . C2H5 has not yet been detected 
in the CSE of CRL 618, however, C2H4 has been detected and the fractional abundance has 
been measured as ~10~ 4 . The calculated fractional abundance of CH3CHO based on the 
formation model is ~5xl0~ 6 . The calculated upper limit based on our observations is below 
this value, X C h 3 cho 10 6 . The implication for the formation and detection of CH3CHO is 
that if the abundance of C2H5 is similar to C2H4 gas phase reactions like those listed above 
are not important to the formation of this molecule. 

One obvious difference between PPNe and molecular cores that may account for the 
difference in grain surface reactions is their different grain properties. In the cold collapse 
phase of molecular cores, molecules are depleted onto ices, where hydrogenation then slowly 
changes the chemical composition of the mantles. The mantles are then evaporated during 
the hot core phase and the subsequent release of these mantle species into the gas phase can 
then drive the complex chemistry observed in hot molecular cores. The envelopes of AGB 
stars and PPNe on the other hand are the actual sites of the formation of carbon dust via 
the formation of PAHs. Hence grain properties are very different from star forming cores 
and the cold collapse phase is completely lacking. Therefore, if grain chemistry and the 
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subsequent processing of the grain mantle products in a hot gas phase is the main route to 
produce large organic biomolecules in a large abundance, they cannot be expected to exist 
in PPNe. For further reviews of the chemistry leading to organic molecules see Tielens & 
Chamley (1997) and Ehrenfreund & Chamley (2000). 

Indeed, while a complex chemisty is taking place toward PPNe, including CRL 618, it 
appears that the physical conditions and chemistry may not be efficient enough to produce 
and detect large organic biomolecules. 


5. Summary 

We have presented the results of our survey toward the protoplanetary nebula CRL 
618 for several large, highly saturated, organic molecules of biological importance includ- 
ing acetaldehyde (CH 3 CHO), acetic acid (CH 3 COOH), dimethyl ether (CH 3 OCH 3 ), ethanol 
(CH 3 CH 2 OH), formic acid (HCOOH) and methyl formate (HCOOCH 3 ); large carbon chain 
molecules including methyl cyanide (CH 3 CN), methylcyanoacetylene (CH 3 C 3 N), cyanoacety- 
lene (HC 3 N), cyanodiacetylene (HC 5 N) and C 6 H; and finally smaller molecules including 
34 SO, S0 2 , OC^S and MgNC. No biologically important organic molecules were detected 
however, based on the 1 a rms noise level in the spectral windows containing our large 
molecular species, we were able to place upper limits on the column density and fractional 
adundances relative to H 2 (Table 4). We made the first interferometric detection of CH 3 CN 
toward CRL 618. Using a rotational temperature diagram of the J — 6 — 5 rotational tran- 
sitions, we found a temperature of 88(44) K and from the temperature and distribution of 
CH 3 CN toward CRL 618, we find CH 3 CN is formed in the outer envelope surrounding the 
UC HII region. We also made the first interferometric detection of vibrationally excited 
HC 5 N toward this source as well as detected transitions of the 1 / 7=3 and 4 states of HC 3 N. 
The P-Cygni line profile and corresponding channel maps of vibrationally excited HC 5 N sup- 
ports its distribution in the extended circumstellar envelope expanding away from the central 
star. The detection of vibrationally excited HC 3 N confirmed the temperature structure and 
column density of HC 3 N in the extended envelope found by Wyrowski et al. (2003). Finally, 
our observations clearly indicate that CRL 618 is a good source of large carbon chain species 
but is a very poor source to detect or produce organic species of biological importance. 


We acknowledge support from the Laboratory for Astronomical Imaging at the Univer- 
sity of Illinois, and NSF AST 99-81363 and AST 02-28953. 
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Fig. 1. — Spectra of the four transitions of vibrationally excited HC 5 N located in the pass- 
bands of our survey, (a) First detection the HC5N J=40-39, (0,0,0, 1) 1+ transition toward 
CRL 618 seen in absorption and emission. The emission peaks near a v LS r of -24 km s -1 . 
The spectral line label in the upper left of each panel is the rest frequency of each HC 5 N 
transition for a v LS r of -27.5 km s -1 . The dashed line is centered on the vlsr of CRL 618 
(-27.5 km s -1 ). Near 40 km s -1 , we begin to detect the 34 SO 2(3)-l(2) transition in emission. 
We also label the location of the CgH 2 n 3 / 2 J=77/2-75/2 f transition. The 1 a rms noise 
level is ~40 mJy beam -1 (indicated by the vertical bar at the left of each spectrum), (b) 
Passband containing the HC5N J=40-39, (0,0,0,3) 1- transition toward CRL 618. No spec- 
tral feature is detected beyond the 1 o rms noise level, (c) Passband containing the HC5N 
J=40-39, (0,0,0,3) 1+ transition toward CRL 618. No spectral feature is detected beyond the 
1 a rms noise level, (d) Passband containing the HC5N J=40-39, (0,0, 0,4) transition toward 
CRL 618. We also label the location of the MgNC N = 9 — 8 transition. No spectral feature 
is detected beyond the 1 a rms noise level, however the MgNC transition is beginning to be 
resolved out of the noise (§4.4). 

Fig. 2. — Spectra of the detected transitions of CH 3 CN and HC3N found from our survey 
toward CRL 618. (a) First detection of the J = 6 — 5 rotational transitions of CH3CN seen 
entirely in emission. Between the K — 2 and K = 3 lines of CH3CN lies the J — 12 — 11 
(0,1,0,3) 1+ vibrational transition of HC 3 N seen in absorption. The measured vlsr of the 
absorption trough of HC3N is -27.5 km s -1 . The spectral line label in the upper left of the 
panel is the rest frequency of the HC 3 N transition. The 1 a rms noise level is ~40 mJy 
beam -1 (indicated by the vertical bar at the left of each spectrum), (b) Spectrum of 3 
vibrationally excited transitions of HC 3 N. Each of these transitions lie at a vlsr of -27.5 km 
s -1 . The spectral line labels located in the upper left of the panel corresponds to the rest 
frequency of the HC3N J = 12 — 11 (0,0, 0,4)° transition. 

Fig. 3. — HC5N velocity-channel maps toward CRL 618. The data displayed are a combina- 
tion of the B and C array data. The LSR velocities of the absorption and emission features 
are indicated in km s -1 units at the top right of each panel. The contour level ranges are 
-0.24, -0.20, -0.16, -0.12, -0.08, 0.08, 0.12, 0.16, 0.20, 0.24 and 0.28 Jy beam -1 . The syn- 
thesized beam of 5."8x5."l is shown at the lower left of each panel. The last panel shows 
the 3 mm continuum emission toward CRL 618 in the BIMA B array The synthesized beam 
of 3. "0x1. "9 is shown at the lower left of this panel. The continuum emission peak at an 
intensity of 1.86 Jy beam -1 at a frequency near 108 GHz. 

Fig. 4. — HC 3 N v 7 =4 velocity-channel maps toward CRL 618. The data displayed are a 
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combination of the B and C array data. The panel labels and beamsizes are the sam p as in 
Figure (2). The contour level ranges are -0.36, -0.32, -0.28, -0.24, -0.20, -0.16, -0.12, -0.08, 
0.08, and 0.12 Jy beam -1 . 

Fig. 5. — Rotational temperature diagram for the observed CH 3 CN transitions toward CRL 
618. Ordinate: In (N u /g u ). Abscissa: E U (K). The least- squares fit is shown as a solid line. 
The resulting rotational temperature is 88(44) K in a beam size of 5."8x5."l. 
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Table 1. Observational Parameters 


Observation Date 

Gain Calibrator 

Flux Calibrator 

Frequency 

Array 

Beam 

(i) 

(2) 

(3) 

(GHz) 

(4) 

(5) 

(arcsecxarcsec) 

(6) 

2002 Oct 11 

0359+509 

Saturn 

1X1.5 

c 

7.18x6.59 

2002 Oct 19 

0359+509 

Saturn 

111.5 

c 

6.97x6.50 

2002 Nov 04 

0359+509 

Jupiter 

110.3 

c 

6.91x6.43 

2002 Nov 18 

0359+509 

Saturn 

111.5 

c 

7.46x6.66 

2003 Jan 28 

0359+509 

Jupiter 

110.3 

B 

2.90x1.84 

2003 Jan 31 

0359+509 

Jupiter 

111.5 

B 

3.72x1.90 

2003 Fteb 01 

0359+509 

Jupiter 

111.5 

B 

2.75x1.96 

2003 Ffeb 03 

03-59+509 

Jupiter 

111.5 

B 

3.48x2.04 



Table 2. 

Molecular Line Parameters 



Species 

Transition 

Frequency* 

(MHz) 

< Sijf > 
(Debye 2 ) 

E» 

(K) 

Qr 

C 6 H b 

2 n 3/2 J=77/2-75/2 

f 106,762.470(90) 

1340.2 

101.2 

89.8T 

M SO c 

2(3)-l(2) 

106,743.374(18) 

3.6 

21.0 

2.9T 

OC 34 S d 

9-8 

106,787.380(80) 

4.6 

25.7 

3.4T 

so 2 b 

273,2S'264,22 

107,060.209(3) 

8.3 

370.5 

1.1T 3 / 2 

MgNC b 

9-8 

107,399.800(400) 

245.2 

25.9 

7.0T 

ch 3 cn« 

6(4)-5(4) 

110,349.760(200) 

51.0 

132.3 

1.0T 3 / 2 


6(3)-5(3) 

110,364.470(100) 

68.9 

82.6 



6(2)-5(2) 

110,374.986(120) 

81.6 

47.0 



6(1)-5(X) 

110,381.376(120) 

89.3 

25.7 



6(0)-5(0) 

110,383.494(120) 

91.9 

18.6 


(CH 3 ) 2 O f 

9l,8-82,7 AA 

106,775.679(80) 

3.66 

43.4 

91.6T 3 / 2 


9i , 8 - 82,7 EE 

106,777.371(80) 

3.66 

43.4 



9i,8-82,7 AE 

106,779.058(80) 

3.66 

43.4 



91 , 8 - 82,7 EA 

106,779.066(80) 

3.66 

43.4 


CH 3 CHO b ’ 9 

3i,2-2o,2 E 

107,887.477(1598) 

0.02 

8.07 

2 3t 3 /2 

CH 3 COOH h 

10*,io-9*,9 E 

111,507.270(40) 

54.8 

30.5 

14.1T 3 / 2 


10*, io-9*, 9 A 

111,548.533(40) 

54.8 

30.5 


CH 3 C 3 N‘ 

27(6)- 26(6) 

111,503.905(3) 

579.1 

345.3 

8.7T 3 / 2 


27(5)- 26(5) 

111,515.692(3) 

588.3 

262.8 



27(1)-26(1) 

111,541.424(4) 

608.4 

82.7 



27(0)-26(0) 

111,542.497(4) 

609.2 

75.2 


ch 3 ch 2 oh* 

127,6-136,7 

111,510.340(31) 

1.7 

127.4 

6.5T 3 / 2 


127,5-136,6 

111,510.550(31) 

1.7 

127.4 



82 , 7 - 80,8 

111,738.020(19) 

0.1 

35.3 


HCOOH k 

5o,5-4o,4 

111,746.780(10) 

10.0 

16.2 

1.7T 3 / 2 

HCOOCHs 1 

97 , 2 - 87,1 E 

110,525.598(32) 

59.2 

9.6 

12.4T 3 / 2 


97 , 3 - 87,2 A 

110,535.182(28) 

59.2 

9.6 



97 , 2 - 87,1 A 

110,535.184(28) 

59.2 

9.6 



97 , 3 - 87,2 E 

110,535.955(26) 

59.2 

9.6 



72,6-6i,5 E 

110,550.217(28) 

19.0 

1.4 
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Table 2 — Continued 


Species 

Transition 

Frequency* 

(MHz) 

< > 
(Debye 2 ) 

Ett 

(K) 

Qr 


72,6-6i,5 A 

110,560.053(30) 

19.0 

1.4 



9i,s-8i,7 E 

111,674.070(13) 

22.8 

28.9 



3i , 8-81,7 A 

111,682.224(13) 

22.8 

28.9 



10i,io-9o,9 E 

111,733.936(14) 

3.8 

30.3 



10i,io-9o,9 A 

111,735.329(15) 

3.8 

30.3 


HC 3 N m 

J=12-ll (0,0, 0,4)° 

110,543.700(130) 

169.4 

1314.4 

13.7T 


J=12-ll (0,0,0,4) 2_ 

110,549.013(1594) 

169.3 

1314.4 



J=12-ll (0,0, 0,4) 4 

110,554.377(42) 

169.3 

1314.4 



J=12-ll (0,1,0,3) 1+ 

110,366.474(337) 

169.0 

994.4 


HC 5 N n 

J=40-39 (0, 0,0,1)’+ 

106,767.4842(44) 

749.96 

259.4 

15.7T 


J=40-39 (0, 0,0,3) 1_ 

107,068.1920(300) 

750.0 

568.2 



J=40-39 (0,0, 0,3)’+ 

107,254.3370(300) 

750.0 

568.4 



J=40-39 (0,0,0,4) 4 - 

107,394.3161(50) 

750.0 

726.2 



J=40-39 (0,0,0,4) 4 + 

107,394.3435(50) 

750.0 

726.2 



a Errors are 2 a 
b Pickett et al. (1998) 
c Tiemann (1974) 
d Maki (1974) 

'Boucher (1980) 
f Groner et al. (1998) 

SBauder (1976) 
h Illuy8hin et al. (2001) 

'Muller et al. (2001) 

JLovas (1982) 
k Willemot et al. (1980) 

1 Oesterling et al. (1999) 

“Wyrowski et al. (1999), Yamada k C reswell (1986) 
"Yamada et al. (2004), Muller et al. (2001) 
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Table 3. Line Detections 


Frequency 

(MHz) 

Species 

TYansition 

vlsr 
(km s _1 ) 

/ Aldv 

(Jy beam -1 km s -1 ) 

106,767.4842(44) 

hc 5 n 

J=40-39 (0,0,0,1) 1+ 

-24 

4.74(1.9) 




-34 

-4.15(1.7) 

106,743.374 

34 SO 

2(3)-l(2) 


& 

110,349.760(200) 

CH 3 CN 

6(4)-5(4) 

-27.5 

0.58(42) 

110,364.470(100) 

CH 3 CN 

6(3)-5(3) 

-27.5 

2.77(18) 

110,366.742(337) 

hc 3 n 

12-11 (0,1,0,3) 1+ 

-27.5 

-2.69(30) 

110,374.986(120) 

CH 3 CN 

6(2)-5(2) 

-27.5 

1.85(17) 

110,381.376(120) 

ch 3 cn 

6(1)-5(1) 


b 

110,383.494(120) 

CH S CN 

6(0)-5{0) 

-27.5 

7.63(24) b 

110,543.508(130) 

hc 3 n 

12-11 (0, 0,0,4)° 

-27.5 

-1.56(3) 

110,550.328(1594) 

hc 3 n 

12-11 (0,0, 0,4) 2- 

-27.5 

-1.38(3) 

110,554.742(42) 

hc 3 n 

12-11 (0,0, 0,4 ) 4 

-27.5 

1.92(3) c 


a Because we did not see the line peak, no attempt was made to fit this transition. 

b Because the K — 0 & 1 components could not be fit with separate Gaussians, a one Gaussian 
fit was applied to both lines. 

c This line may be blended with the 154,2-145,2 transition of C-C 3 H 2 at 110.558(9) GHz (Pickett 
et al. 1998). 
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Table 4. Column Densities 


Species 

Frequency 

Nrot 

X(H 2 ) 


(MHz) 

(cm -2 ) 



Detections 

hc 5 n* 

106,767.484 

9.9 xlO 14 

7.6x10“® 

hc 3 n 

b 

4.7(6) xlO 17 

3.6(5) xlO -6 

ch 3 cn 

b 

3.2(3) xlO 14 

2.4(3) x lO" 9 

Upper Limits 

C 6 H 

106,762.470 

2.5X10 13 

1.9 xlO" 10 

OC M S 

106,787.380 

1.2xl0 14 

9.1 xlO" 10 

so 2 

107,060,209 

l.OxlO 16 

7.8xl0“ 8 

MgNC 

107,399.800 

4.4 xlO 12 

3.4xl0 -11 

HCOOCHa 

110,550.217 

1.1 xlO 16 

8.5x10-® 

(CH 3 ) 2 o 

106,777.371 

l.OxlO 17 

7.9xl0“ 7 

ch 3 cho 

107,887.477 

1.4xl0 17 

1.1 xl0“ 6 

ch 3 cooh 

111,548.533 

3.5 xlO 1 4 

2.7x10-9 

ch 3 c 3 n 

111,542,497 

3.3xl0 13 

2.5x10- 10 

ch 3 ch 2 oh 

111,738.020 

9.4xl0 16 

7.3 xlO" 7 

HCOOH 

111,746.780 

2.0X10 14 

1.5X10- 9 


& Column density based on the emission line detection. (§4.3) 

b All lines were used to calculate the total column density of 
this molecule. 


